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FFECTIVE communication is essential to technical prog- 
ress. This is particularly applicable to the ever changing 
field of electronics. The rapid advances tax the capacities of 
all of us to keep pace. Raytheon’s problem in this connection 
is typical — we must keep each other informed of our tech- 
nical activities. In recognition of this need, ELECTRONIC 
PROGRESS was established five years ago as a link in the 
communications chain between Raytheon engineers and scien- 
tists. The growth and acceptance of the magazine attest to the 
importance of its role. 

In assuming the position of Vice President, Engineering and 
Research, I should like to endorse most strongly the view that 
effective communication is vital to our success. A large system 
development program, for example, requires the combined ef- 
forts of many people, from separate organizations, of diverse 
skills and interests. It is imperative that our many working 
groups draw upon each other’s resources, to reap the benefit of 
prior experience, to prevent unnecessary duplication, and to 
lessen the possibility of seeking outside a skill that already 
exists within. For these reasons I urge you to communicate with 
each other, to exchange reports, to organize seminars, to use 
ELECTRONIC PROGRESS as your journal. 

And while we are strengthening our links within, we should 
not neglect our ties to the scientific community outside. It is 
important that we frequently publish in the professional and 
trade press, and regularly report our findings at the meetings 
of the societies. This is both a responsibility and a privilege 
since there is no better way to professional advancement. 

Your cooperation in these efforts will produce a greater 
understanding of what Raytheon is and does — not only among 
ourselves but among our colleagues in the industry at large. 


Martin ScHILLING 
Vice President, 
October 19, 1960 Engineering and Research 












NTRODUCTION 


The theme of this issue, ‘‘ Mathematics in Indus- 
try,’’ is an ambitious one. It would be impossi- 
ble in one issue to cover the broad impact of 
mathematics on industry even limiting our in- 
terest to electronics. Our purpose, however, is 
a much narrower one — to acquaint the engineer 
with some of the mathematical tools available to 
him. Computers and statistics are presented as 
prime examples. In future issues some of the more ad- 
vanced developments in mathematical approaches to 
engineering problems — operations research, informa- 
tion theory, learning machines — will be discussed. 

Some comparison of the analog and digital com- 
puters and their particular applications would perhaps 
be helpful. 

The analog computer essentially presents a visual 
working model of the mathematical function being 
solved. Physical quantities (voltages in electronic 
analogs, shaft positions and light intensities, for ex- 
ample, in other forms of analogs) are forced to obey 
mathematical relations comparable to those describing 
the physical system being analyzed. The analog com- 
puter operates in ‘‘real’’ time; i.e., the analog process 
oceurs during the same time as does the variation of 
the parameters it represents. The output of the ana- 
log is thus a continuous display of the variation of the 
parameters, and is conveniently recorded, either per- 
manently or temporarily, on strip charts, plotters, 
oscillographs, ete. In addition, in most analog com- 
puters, each step or operation of computation is per- 
formed by a separate unit, and all units operate 
simultaneously. This parallel operation makes the 
solution available almost immediately. 


In the electronic digital computer, on the other 
hand, data is not processed in a form analogous to the 
input information, but rather in the discrete form of 
a series or set of characters (language form). Because 
of the basic operation of its internal computing ele- 
ments (flip-flops, counters, magnetic core storage ele- 
ments, ete.), which can represent only two states, on 
or off, the digital computer operates in the binary 
number system. The computer consecutively performs 
arithmetical operations in a predetermined sequence 
rather than in the paralleled operation of the analog. 


In contrast to the analog computer, which can not 
perform operations that are not readily representable 
in analog form (some nonlinear functions), the digital 
computer can perform any type of mathematical opera- 
tion. In terms of accuracy, the analog is limited by the 
precision of its components ; the accuracy attainable in 
most practical analog installations is limited to 0.01 
per cent or less. With the digital computer any degree 
of accuracy can be obtained and maintained if a suffi- 
cient number of characters are used in the language 
code of the system. In straight binary code, for ex- 
ample, accuracies of the order of 0.001 per cent can 





be obtained if ten decision elements are used, or 0.0005 
per cent if eleven decision elements are used. Hence, 
the degree of accuracy attainable depends directly on 
the capacity of the digital installation and not on the 
precision of the components. 

Due to the construction of the machines the formu- 
lation of a program for each is necessarily different; 
the original preparation of a problem, however, is gen- 
erally the same for both. The mathematical relation- 
ships, boundary and initial conditions, numerical val- 
ues of constants, range of parameter variations, and 
estimated results must be determined before the actual 
programming process is begun. 

After the initial analysis, the programming proce- 
dure for the analog involves the scaling of equations 
so that the maximum excursions of the variables do 
not exceed machine ratings. The next stage is the 
formulation of the computer block diagram for each 
equation. The block diagram specifies all computing 
elements to be used, necessary interconnections, set- 
tings of initial conditions and constants, ete. The last 
step is the actual wiring of the anaiog patch board in 
accordance with the block diagram. 

In programming a digital computer a set of detailed 
instructions must be formulated in terms similar to 
the language of the computer. Since the computer 
will perform only one operation at a time, these in- 
structions must be written in detail, one after another, 
to provide the procedure by which the machine will 
perform the solution of the problem. 

The instructions are typically written first in flow 
diagram form, in which a graphical picture of the 
course of the problem’s solution is given. When there 
are no decisions to be made by the computer control, 
the flow diagram will be a continuous line interrupted 
only by boxes indicating the sequence of operations in- 
volved. When iterative calculations or repetitive data 
processing sequences are required loops will appear in 
the flow diagram. Depending on the requirements of 
the specific problem the diagram may be very detailed, 
or may show only a major outline of the problem, with 
dozens of machine instructions implied in each box. 
The process of translating the instructions given in 
the flow diagram boxes into sets of specific machine- 
language instructions is called coding. Coding may be 
done either by hand, or by means of one of the auto- 
matic coding systems described below. 

The choice between an analog or a digital computer, 
or some other means of computation, will depend on 
the requirements of the problem to be solved. How 
complex is the problem? Is it adaptable to the analog 
computer? Are the programming time and expense for 
the digital computer justified? Is the computation 
repetitive, in which case the digital computer has a 
special advantage? No set rule can be formulated as 
to when to use a desk calculator, or automatic equip- 
ment, or which type. It is hoped that the following 
articles will make the engineer’s task in making a 
choice somewhat simpler. 
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ANALOG 
COMPUTATION 


A S early as 1916, electric problem boards were used 
by the electric power companies to predict the opera- 
tion of their distribution systems when a fault (short 
or open circuit) occurred on one of the lines. These 
boards were lumped-resistance, inductance models of 
the distribution system being studied with batteries 
or small AC generators representing the generating 
stations. Ammeters and voltmeters were used to 
measure the effects of a fault anywhere on the lines. 
By 1930 these plotting boards had been refined to 
the point where line voltage stability with dynamic 
load variation could be studied.’ ? 

Alternating current plotting boards are still being 
used by the power industry. Modern machines are 
electronic in nature and operate at frequencies as high 
as 10 ke. Three-phase systems with variable phase 
control of the generators and loads can be analyzed. 
These plotting boards are special purpose analog ma- 
chines and can not be used for solving generalized 
differential equations.* 

Between 1927 and 1931 Vannevar Bush developed 
the first large scale machine for the solution of dif- 
ferential equations. This machine, the first differen- 
tial analyzer, used wheel and disc integrators and later, 
servomechanisms. Although very slow by modern 
standards, it was capable of five-figure precision. 
Later models of this machine were used until the ad- 
vent of the modern analog computer during and after 
World War IT. 

The Mechanical Transients Analyzer is a type of 
analog computer developed by Westinghouse during 
World War II. In this simulator a mechanical sys- 
tem is studied by measuring the response of an elec- 
trical system that has an analogous set of differential 
equations describing its operation. Although the ma- 
chine was designed to analyze mechanical systems, it 
can handle any set of differential equations that has 
an analogous electrical system.* 

The modern electronic analog machine was made 
possible by the development of operational amplifiers 
during World War II. The operational amplifier is 
usually credited to Lovell of Bell Telephone Labs, 
although such a device was used as early as 1938 by 
Philbrick in computing circuits for the solution of 
servomechanism problems.* Modern machines made 
possible by the operational amplifier are known com- 
mercially by such names as Boeing, Gap/R (George 
A. Philbrick Researches, Inc.), GEDA (Goodyear), 
EASE (Beckman), and REAC (Reeves). 


-_— 


*Philbrick’s work was never published. 
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ALLAN R. THAYER is a Senior En- 
gineer in the Ordnance Servo Department, 
Surface Radar and Navigation Operations. 
His analog equipment experience includes 
several years at Sperry Gyroscope Co. where 
he analyzed the kill probability of an air-to- 
air guided missile using a REAC. Since join- 
ing Raytheon in 1954, he has worked on 
servo design for the Hawk BCC, the SESCO 
range converter (a digital-to-analog con- 
verter), and the AQS-8 transducer drive. 
Mr. Thayer received his B.S. and M.S. de- 
grees from Worcester Polytechnic Institute 
in 1952 and 1954, respectively. 


JOHN P. WARD is a Design Group 
Leader in the Ordnance Servo Department, 
Surface Radar and Navigation Operations. 
Mr. Ward designed Nautilus reactor con- 
trols at Westinghouse for two years before 
joining Raytheon in 1954. At Raytheon he 
has designed servo systems for Hawk, SPG- 
5i, AQS-8, and SESCO. He also worked 
with digital and analog equipment during 
this period. 

Mr. Ward received his B.S. degree in EE 
from MIT in 1952. He has published two 
papers in the Transactions of the AIEE. 





Electronic Analyzer Building Blocks 


A block diagram of the generalized operational 
amplifier, a high-gain DC amplifier designed for low 
drift, is shown in Figure 1. Z;,, the feedback im- 
pedance, determines the type of operation performed, 
while Z,, Zo, . . . Z, are all resistive elements. 

The three most common types of operational ampli- 
fiers are the summer, the inverter, and the integrator. 
Figure 2 shows the block diagram representation of 
each of these three operational amplifiers. In the sum- 
mer, the feedback impedance is a resistor and the 
input impedance is composed of resistors of various 
values. The magnitude of the ratios of Rr/R,, Rr/Reo, 
ete. determines the gain associated with each particu- 
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Figure 1. Generalized Operational Amplifier 
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Figure 2. Block Diagrams for Summer, Inverter, and 
Integrator 


lar input. The inverter, which is usually used as a 
sign changer, is basically a summer with a gain of 
one, since its feedback resistance is always equal to 
each input resistance. In the integrator, the feedback 
impedance is a capacitor and the input resistors may 
be varied in magnitude to obtain various gains, as 
with the summer. 

A differentiator could be made using capacitive in- 
put impedances and resistive feedback. Such an 
operational amplifier would be very sensitive to noise, 
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however, because the gain of a differentiator increases 
with frequency. As an example, a differentiator for 
computer such as a REAC (Reeves Electronic Analog 
Computer) would be very difficult to shield adequately 
from 60-cps noise because the normal problem fre. 
quencies are between 1 and 10 eps. Thus for normal 
operation, 60-cps noise going through a differentiator 
would be effectively multiplied by a factor of from § 
to 60. There is, fortunately, no need for differentia. 
tion in analog computer work since differential equa- 
tions can be solved with the use of integrators instead 
of differentiators. 

Very complicated transfer functions, those involy- 
ing filters for servomechanisms, for example, can be 
represented by one operational amplifier with ap. 
propriate input and feedback impedances. Transfer 
functions involving ratios of eighth-order differential 
equations in ‘‘s,’’ the Laplace transform variable, are 
available from one DACL (Dynamic Analysis and 
Control Laboratory) operational amplifier. 

With the devices discussed thus far only linear sys- 
tems may be solved. If a multiplying unit is available, 
however, many nonlinear systems may be studied. 
Two of the several types of multipliers available for 
analog computers are shown in Figure 3. 

Figure 3a is the quarter-square multiplier. A suit- 
able squaring unit commonly used for this multiplier 
is a triode tube operating near cut-off. Since the band- 
width of the quarter-square multiplier is limited by 
the associated summing amplifiers and not the squar- 
ing unit, a large bandwidth is available. The disad- 
vantage of the quarter-square multiplier is that if the. 
input signal varies over even a moderately large 
range of magnitudes, the output departs from the true | 
square. For large input signals the output of the 
squarer becomes a linear function of the input, and 
the device ceases to be accurate. Great care must be 
taken in controlling signal amplitudes in a problem 
handled by a quarter-square multiplier. Signal levels 
must be such that the triode operates near cut-off and 
the operational amplifiers operate at voltage levels and 
in voltage ranges for which they are designed. This 
control of voltage levels and ranges in a given prob- 
lem is often referred to as ‘‘amplitude scaling.’’ 
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3a. Quarter-Square Multiplier 
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Figure 3b shows the block diagram of a servo mul- 


tiplier. To multiply two variables in this device a 


servo-driven potentiometer is excited by one variable, 
and the shaft is positioned by the other variable. The 
voltage read at the wiper is proportional to the prod- 
uct of the two variables. In direct contrast to the 
quarter-square multiplier, the servo multiplier can 
operate over a wide range of inputs, but has a limited 
bandwidth. The frequency limitation of the servo 
multiplier in the REAC installation at the Wayland 
Laboratory is approximately 1 cps at 75 volts peak 
input amplitude; at 10 volts peak amplitude and 1.5 
eps the servo error is less than 0.5 per cent. 

In the REAC installation, six servo-driven poten- 
tiometers are ganged on one shaft. One potentiometer 
is used by the servo, and is not available for comput- 
ing. Three of the other five are available for multiply- 
ing. The remaining two potentiometers have taps, and 
are usually used for function generation by appropri- 
ate loading of the taps. The tapped potentiometers 
also may be used for multiplying, if necessary. 

Division of two variables also is possible with servo- 
driven potentiometers by appropriate connection pro- 
vided the divisor variable is always non-zero, since no 
machine can handle division by zero. Figure 3c shows 
the block diagram for a servo divider in which Z = 


100 +, where the 100 is the servo seale factor. The 


computer solves the relation : 


[x+in [—A] =Z, or Z= A 
00 71 
— 100X 
ee ,for $Y» 1. 


For values of Y above 1.0 volt this division is accu- 
rate on the REAC, but for values of Y lower than 
1.0 volt, accuracy is degraded. 

DC resolver servos utilizing sine-cosine poten- 
tiometers may be employed to resolve or generate 
functions. Two orthogonal vectors can be easily con- 
verted into R,® polar form. The resolver servo is 
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3b. Servo Multiplier 


capable of following slow changes in the input vari- 
ables with a dynamic error of less than 0.5 per cent 
of the applied voltages at 0.8 cps with 20 volt peak- 
to-peak inputs. 

Another important nonlinear generating element is 
the diode function generator (DFG). When using the 
DFG any arbitrary function with respect to an inde- 
pendent variable is approximated by a series of 
straight-line segments. Each segment has a break- 
point and slope. The transfer function for each seg- 
ment is: 


xX,=m (Xi—p), 


where X, = output, X; = input to DFG, 


m = slope, and p = breakpoint. 

In the REAC it is possible to have slopes of +12, 
and breakpoints that can vary from —100 volts to 
+100 volts. The slope of successive segments is equal 
to the algebraic sum of the slopes of each segment that 
has breakpoints less than the DFG input. The tran- 
sitions of the diodes from non-conduction to conduc- 
tion are gradual, making the breakpoint regions 
slightly rounded. 

The free-diode limiter is also an important non- 
linear generator, and can generate functions which 
simulate backlash, dead-zone, and saturation. 

Operational relays can be used to switch variables, 
factors, or conditions. With relays a problem may be 
held static to measure parameters, or stopped at pre- 
set boundary conditions. 


Analog Output Data Presentation 


Output data is presented by the computer in vari- 
ous forms; the more common ones are cathode-ray 
oscilloscope graphs, multiple-channel strip recordings 
and two-axis plots. The oscilloscope is used with 
high-speed analog machines, while strip recorders 
(rectangular coordinate devices) are used with slower 
analog machines. The two-axis plotter has particular 


utility in presenting cross plots such as aus vs X. 
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Figure 3. Quarter-Square Multiplier and Servo Multiplier and Divider 
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Examples of Analog Computation 


The most important use of the analog computer in 
pure and applied mathematics is the solution of dif- 
ferential equations. The order and type of the equa- 
tion that can be solved is limited only by the capacity 
of the analog installation. Higher-order differential 
equations which are very tedious or even impossible to 
solve analytically can be solved with the computer. 


To illustrate the method of a problem solution con- 
sider the equation : 





d2x(t) +10 dx(t) 


- Sy + 1000x(t) = 35A, 


with initial conditions: 


dx(t) 


dt t= . ile 


and x(t) = —X,. 


t= @ 


The given equation may be transformed as follows: 
s?x + 10sx + 1000x = 35A, 


where s = Laplace operator. The natural frequency 
of this equation is: 


®y = V 1000 = 31.6 rad/sec. 


Because the computer is limited to 10 to 20 eps in its 
frequency range, it is often necessary to effect a time 
scale change between problem time and computer time. 
For example, this problem will be solved with a time 
scale reduction of 10, thereby reducing the frequency 


on the computer to 3.16 rad/see. The conversion equa- 
tions are: 


T=10t T = Computer independent time 
variable 

10SX = sx S = Computer transform vari- 
able 

100S82X =s?x X = Computer dependent vari- 
able. 


The machine equation (using a 1:1 analog for vari- 
ables) becomes: 


100S2X + 100SX + 1000X = 35A 
100S2X = 35A — 100SX — 1000X 
$2X = 0.35A — SX — 10X. 


To arrive at the computer circuit it is assumed that 
the highest-order derivative of the desired function is 
known. This signal is fed through the required num- 
ber of integrators to arrive at the desired output. With 
the output and all derivatives available, appropriate 
feedback to the first integrator satisfies the equation 
for the highest-order derivative. Figure 4 shows the 
block diagram for the solution to the above equation. 
Notice that the dependent variable and its derivatives 
ean all be monitored simultaneously vs the independ- 
ent variable. 

An analog computer is useful in solving equations 
of many forms. Some partial differential equations 
have been solved by converting the problem to one of 
ordinary equations. Parabolic and hyperbolic partial 
differential equations of the second order are solved 
by a difference method. Starting with known function 
values along a characteristic curve, a new characteris- 
tie curve is obtained, and the process continues. 

The real and imaginary parts of a polynomial 


P,(z) = a + ayz +... + agz", where z = re® also 7 


have been evaluated on analog equipment. The roots 
of polynomials were found by making a systematic 
search for values of r and § making P,(z) = 0. 

Other functions performed by the analog computer 
include the evaluation of multidimensional integrals 
and the generation of a vector with N-dimensional 
normal distribution. 

An interesting application of an analog machine is 
the evaluation of the Fourier integral: 


F(@) = fr) edt. 


The magnitude of this function as a function of fre- 
quency, ©, is often called the frequency spectrum of 
f(t). The function f(t) represents the waveform and 
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is assumed real, while F(@) is generally complex. The 
converse problem, determining the waveform f(t) cor-" 
responding to a given frequency spectrum, also arises. 
f(t) can be evaluated by means of the inverse Fourier 
transform : 


f(t) = = i ” F(a) edo, 


A solution has been outlined which consists of de- 
veloping and solving a differential equation whose 
analog solution is or leads to the required function.5 
This differential equation is derived from the Fourier 
integral and the Fourier series for an arbitrary func- 
tion of time and, for time greater than zero and initial 
conditions equal to zero, takes the form: 


Srtt) +02 f yt) dt = f(t). 


The analog representation is shown in Figure 5. 


The electronic analog yields the real part of the 
solution, y(t). Therefore, for any numerical value of 
@, the electronic analog solution gives the correspond- 
ing magnitude and phase of the direct Fourier trans- 
form of the arbitrary function f(t). 

Physically, the equation describes the response of a 
simple oscillatory system due to an arbitrary excita- 
tion defined for t>0. The steady-state energy of the 
system, which oscillates at an angular frequency , 
is absorbed from the excitation. The oscillation ampli- 
tude is a measure of the relative strength of the 
Fourier components of f(t) at the corresponding fre- 
quency. 

Complete radar RF and video portions have been 
simulated on analog computers.* Simulation of radar 
systems has aided in the analysis of radars and missile 
guidance systems and in evaluation of radar vul- 
nerability to countermeasures. Such simulation has 
enabled the system: designer to evaluate range resolu- 
tion, velocity resolution, sensitivity, noise perform- 
ance, and vulnerability to jamming. 
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Figure 5. Solution of Fourier Integral 
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Figure 6. Radar Simulation: 


In simulating ‘radar systems, particular attention 
must be paid to time sealing. Thus, although the scal- 
ing of RF and IF center frequency is somewhat arbi- 
trary, all transient bandwidths (RF, IF, video, and 
servo) in a radar system must be scaled by the same 
scale factor. The IF center frequency is usually scaled 
to be about one-millionth of the system IF. Typical 
intermediate frequencies on the computer are 10 to 
30 eps. Since the system RF carries no information, 
it may be scaled to an audio frequency that permits a 
convenient IF ard bandwidth. 
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Some examples of radar processes that can be simu- 
lated are shown in Figure 6. Figure 6a shows the 
exact equivalent to a single-tuned circuit. Staggered 
pairs and triples can be simulated by cascading the 
single-tuned circuit. Figure 6b shows the computer 
simulation of the linear and square-law detectors, and 
Figure 6e shows a simulated delay line on which one 
tape cycle is the time delay between line taps. 

Problems in the fields of heat transfer, hydraulics, 
vibration, circuit analysis, flight dynamics, and others 
where differential equations characterize the system 
are usually adaptable to analog computer solution. 
Problems characterized by a pulse waveform such as 
sampled-data or radar systems are also adaptable to 
computer simulation. 


Possible Accuracies 


Analog computer solutions are obtainable to accu- 
racies of 0.1 per cent. When accuracies better than 
1 per cent are required, the cost compared to the cost 
of a digital computer becomes prohibitive. In many 
engineering applications solution accuracies of 5 per 
eent are quite adequate when one realizes that the 
original data to be programmed is no better. 

Raytheon’s REAC installation can provide accu- 
racies of 1 to 5 per cent depending on the size and 
nature of the problem. Although the computer inte- 
grators, summers, and inverters are operable up to 20 
eps, problems that require multiplication, division, and 
resolving must be implemented by servos. The servo 
frequency range is limited to about 1.5 eps where the 
servo error becomes 0.5 per cent of input for a 10-volt 
peak amplitude sinusoidal input. 
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INVENTIONS 


Congratulations to the following persons for their valuable contribution 
to Raytheon’s portfolio of inventions: 


Name Field of Invention 


Raymond R. Nelson. ses+ses+++.. Magnetostriction Filter 
Mason W. Huse............. 


Rectifier Meter Circuit 
Design for Super-Power Tube 
. Oil Well Microwave Tool 
MU Mis PI AIS. on 0 5s. canccssmsscowaan Band for ‘“Taut-Band’’ Meter 
Melvin Prager High Speed Analog Digital 
Clark F. Crocker Computer 








At the present time the Wayland Laboratory com- 
puter facility has a REAC with capacities noted in 
the accompanying table. By October 1, 1960, the 
facility will be increased by the addition of an analog 
machine obtained from the Dynamie Analysis and 
Control Laboratory (DACL) at MIT. The addi- 
tional capacity added to the installation is also noted 
in the table. The DACL setup is a console in which 
modules are easily changed to implement functions 
from a stock of spare modules. 


The DACL and REAC machines will be set up so 
that they may be operated either independently on 
separate problems or in cascade on the same problem. 
The over-all installation will have roughly the capaci- 
ties of a computer rated at 100 operational ampli- 
fiers. The Wayland computer is available for use by 
anyone in Raytheon. For ‘problem analysis or con- 
sultation consult Allan Thayer or John Ward at the 
Wayland Laboratory. 


WAYLAND ANALOG COMPUTER 
FACILITY 


REAC installation 


Operational amplifiers normally used as 
12 Integrators (convertible to summers or in- 
verters) 
12 Summers 
16 Inverters (convertible to summers or inte- 
grators) 
Seale factor potentiometers 
Servos 
2 Multiplying servos 
2 Resolving and multiplying servos 
Function generator channels consisting of 
2 Diode function generator channels 
2 Channels of loaded potentiometer function 
generators on each servo 
Free diodes 


Function relays 


DACL installation 


Operational amplifiers usually used as 
9 Integrators 
10 Summers 
8 Lag networks 
3 Quadratic lag networks 
12 Coefficient networks 


Scale factor potentiometers 


Servos each of which can be used for any of the 
following operations 

Integration 

Resolving 

Multiplication 
Function generators 
Diode clipper networks 
Delay units (0.1 sec to 2 sec.) 
Square root extractors 
Squarers 
Sample and hold units 
Sampling pulse generator 


Function relays (each relay has 4 sets of SPDT 
contracts) 
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Electronic switches 

44-step coefficient changer 
X, Y recorders 

Z(X, Y) function generators 
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Tue continuing revolution in data processing ma- 
chines is being felt by today’s engineer in many ways. 
The financial applications are apparent to him almost 
daily. He buys gasoline for his automobile by signing 
the card that will be punched and returned to him 
with his statement itemized and totaled by machine. 
His driver’s license and professional engineering 
registration are prepared on a punched card, part of 
which he carries in his billfold. Even his paycheck is 
prepared by a tabulating machine from a punched 
eard record. 

He is undoubtedly also aware of contributions made 
to the solution of complicated design problems in his 
professional field by workers employing large scale 
digital computers. If he has not by now already uti- 
lized data processing machines or computers on his 
own engineering problems, it is becoming increasingly 
probable that he will do so in the near future. The 
improved design, general acceptance, and wider avail- 
ability of machines of all capacities makes this true. 

The progress in the design of machines of prodi- 
gious arithmetic speed and gigantic capacity is only 
part of the story, for the improvement of computers 
of moderate capacity and of the machines which proc- 
ess data information from punched cards are making 
these machines also of greater value to the engineer. 

The evolution of the data processing machine con- 
tinues to progress along two principal avenues. One 
is in the hardware itself. The real impact of the 
transistor as a replacement for the vacuum tube in the 
arithmetic and control circuits, and the use of ferrite 
cores to provide immediate access storage are really 
only now being felt in commercially available ma- 
chines. 

The other avenue is in the programming and use of 
the machine. A great deal of effort has been devoted 
to the problem of reducing the difficulty of translating 
a practical engineering problem to one which can be 
effectively solved by the machine. Several years ago 
this effort was typified by the evolution of SOAP, 
the Symbolic Optimum Assembly Program for the 
IBM 650. This program in effect used the 650 as its 
own clerk, reducing for the programmer the burden 
of selecting storage locations within the machine and 
allowing him to think in a more easily remembered 
symbolic language rather than in the ten-digit in- 
structions of the machine itself. 

In recent years the FORTRAN compiler also has 
evolved, in which an even larger part of the burden 
of program preparation has been placed on the ma- 
chine. The FORTRAN compiler effectively translates 
a description of the problem in almost mathematical 
language into a usable, if not optimum, program for 
solution of the problem. The language of FORTRAN 
is surprisingly similar to the mathematical language 
already familiar to the engineer. 
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Somewhat similar developments applied to ma- 
chines of less extensive hardware can be expected to 
lead to their more universal use. The trend in smaller 
machines would naturally be in the direction of mak. 
ing more of the operations controlled by the universal 
language of the program. 

One of the utilization problems of the 650, for ex- 
ample, was that while most of the operations were 
controlled by instructions in the program, the machine 
still employed a wired control panel to translate card 
input into the memory and the memory into card out- 
put. The instruction manual describing the rules for 
wiring this panel was of the same order of complexity 
as the manual describing the machine language. A 
good many scientific users of the 650 have employed 
only ecards punched in the ‘‘utility board’’ format, 
and consequently have not explored the machine’s 
capability for treating a ‘‘problem’’ in which alpha- 
betiec information must be processed. 

Newer machines .have substantially reduced this 
problem of card format and the accompanying board 
layout problem. In some of the recent machines the 
arithmetic units and the storage are entirely alpha- 
numeric; i.e., accept either numbers, letters, or some 
special symbols, with equal ease. The card format 
board is eliminated because the machine accepts the 
entire contents of a card into the memory without re- 
arrangement. The blank columns that have been a 
bugaboo of the 650 are read and accepted as blank 
eolumns. The necessary rearranging, or editing, is 
accomplished in the memory with the same language 
that is used for other operations. Similarly, the out- 
put problem has been moved in the memory with 
instructions to rearrange and place in reserved por- 
tions of the memory the information that is to be 
printed out on the tabulator or that is to be punched 
into an output card. 

Designing sophisticated and efficient card systems 
may involve extensive analysis prior to the collection 
of data to ensure that adequate means are designed 
into the system to accept all of the possible variations 
of the data involved. Fortunately, a system need not 
be optimum, nor sophisticated, nor efficient to be use- 
ful and valuable. Indeed in some eases the system has 
been valuable in that it dramatically pointed out the 
variability or the uniformity of the data to be en- 
coded. And, as in other eases, the process of preparing 
the problem for machine processing has shown the 
very nature of the problem to be different from that 
expected and, consequently, susceptible to solution by 
quite different means. 

Actually, intelligent use of data processing ma- 
chines is becoming less and less difficult with the 
progress of time. The ability to program many types 
of machines is becoming less important, and it is neces- 
sary to become familiar with only the most recent 
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machines — in quite the same way that the high school 










student today learning to drive is not required to- 
learn the gear shift positions of the 1924 Dodge and 
the 1921 Model T Ford before he can master his fa- 
ther’s automatic transmission convertible. And even 
if they were so required, the youth of today would 
approach the challenges of a bygone era with a lot 
more self-assurance and intuition than the members 
of the skeptical public that witnessed their introduc- 
tion. The role of data processing machines is a more 
accepted part of our lives today — and this in itself 
is a catalyst to the learning process. 

The programming of the IBM 1401 system, for 
example, can be learned by an operator of average in- 
telligence in about five days of instruction. Therefore, 
it is certainly reasonable that an engineer who, because 
of his prior mathematical training, has already been 
trained to accept logical rules, should be able to read 
the instruction manual and to write a workable, if not 
efficient, program. 

In some respects the IBM 650 computer and its mod- 
ern replacement, the 1401, are well suited to the 
problem of data processing as opposed to scientific 
ealculations. Both of these machines are basically card 
input and card output machines with moderate storage 
capacity and, in the case of the 650, moderate speed, 
as compared to the giants such as the 704 or Univac. 
However, when the problem requires limited process- 
ing of data conveniently available on cards, and when 
the results are also reasonably convenient on cards, 
these machines can make valuable contributions. 

The simplest type of engineering problem on which 
punched card equipment can be effectively employed 
is the type requiring only sorting and listing. Con- 
sider a deck of cards describing the characteristics of 
available microwave tubes. A single card describing 
a particular tube type might contain the identifica- 
tion number, manufacturer’s name, operating fre- 
quency, maximum power output, type of construction, 
required anode voltage, and other characteristics of 
interest. 

One useful listing of this deck would be in order by 
the identification number. The cards are sorted into 
this order by passing the deck through a sorter as 
many times as there are digits in the largest identifi- 
cation number. The cards are put into order by the 
machine starting with the lowest order digit first. 

This deck is then run through a tabulating machine, 
such as the IBM 407, for which a control panel has 
been wired to specify where each field on the card will 
appear on the printed page. Even the tabulating ma- 
chine is capable of testing and making decisions. It 
reads two cards during each machine cycle — the one 
which is printed on the page and also the next one to 
be printed. It can be wired to compare the identifi- 
cation number, for example, of these two cards and 
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to double-space if they are different and to single- 
space if they are the same. In this case duplicate cards 
would be quickly picked out on the resulting listing 
by the change in spacing. Other controls sense the 
approaching end of the page and skip to the proper 
line on the next page, print a stored heading, and 
continue the list. 

The list by identification number prepared in this 
manner would be useful in determining the charac- 
teristics of a tube for which the identification number 
is known. 

We may also want a list, however, that would help 
us determine the identification number of a tube with 
a given power output rating at a given frequency of 
operation. We may take the deck as it now stands 
and sort this time by the field containing the operating 
frequency, and the tabulator would then produce a 
list with the types progressing from low to high fre- 
quency. In this case, since our previous sort was by 
identification number, several tubes having the same 
frequency would be listed with the lowest identifica- 
tion number first. It may be useful to have tubes of 
the same frequency listed in order by power output 
rating. In this case the deck must be sorted first by 
the power output field and then by the operating fre- 
quency field before listing. 

Sorting by these two fields in the reverse order 
could be used to produce a list of tubes in order by 
power output. This would be a convenient way of 
identifying tubes of more than a given power output, 
regardless of frequency. 

Such listings of tube characteristics are currently 
prepared by the Advisory Group on Electron Tubes 
at their secretariat at New York University for all 
classified and unclassified tubes. Similar lists of un- 
classified tubes and transistors are also available com- 
mercially. 

One of the first considerations in the use of punched 
card machines for data processing is the form in which 
the data is to be punched and stored in the card. The 
fact that the card, as used in IBM systems, is designed 
to contain only 80 digits, or characters, of information 
forms one of the limits on the amount of information 
that can be conveniently correlated. For this reason 
various expedients are employed to expand the effec- 
tive information utilized. The most universal of these 
is the use of codes. 

The listing of tube characteristics described should 
carry the designation of the manufacturer on the card 
representing each type. A full 80 columns could easily 
be utilized to spell out in full the name and address 
of the manufacturer. But for this particular applica- 
tion, a numeric code could be compiled in the simplest 
case by merely assigning consecutive numbers to 
unique manufacturers as they are encountered in pre- 
paring the catalog. At the same time a manufacturers’ 
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code deck could be prepared containing this code and 
the manufacturer’s name and address. The cross-ref- 
ence lists required are easily prepared from this deck 
—a list by code number, and an alphabetical listing 
by manufacturer’s name. 

Some sensible guesses as to the order of magnitude 
of manufacturers to be encountered must be made to 
determine the size of the field to be allotted on the 
card for the manufacturer code. Three digits allow 
for 1,000 different manufacturers, four digits for 
10,000, ete. 

Even the relatively simple problem outlined here 
seems to introduce the kind of generic elassification 
that the machine processing systems tend to encour- 
age. In the case of the manufacturer code, the question 
becomes, ‘‘ What constitutes a manufacturer for the 
code?’’ Are the various divisions or plant locations 
of General Electric Company to be treated as separate 
manufacturers? What differentiation will be made 
for wholly owned subsidiaries? What about the divi- 
sions of companies which manufacture at several 
locations — should each location be identified sepa- 
rately ? 

While these questions seem to be introduced only 
because of the need for reducing the manufacturer to 
a code, the problem is not circumvented by employing 
full names, since without careful analysis, the same 
manufacturer may be found under several names, or 


late the corresponding efficiency for any number of 
such cards. In this case only multiplication and divi. 
sion need be carried out, and the machines offer little 
in convenience unless the quantity of required 
calculations is high, or if this is but one step in a 
much larger handling of data which is to be accom. 
plished by machine. 

A problem that might be typical of a higher use 
of the machine’s capability is the statistical analysis 
of test data. Again these inputs are cards, each of 
which contains the results of tests on a particular 
unit under similar conditions. There may be as many 
as ten significant observable characteristics per card. 
A relatively simple program can be prepared for an 
IBM 650 to read a deck of such cards, and then ecal- 
culate and punch out the statistics of each of the ten 
characteristics — number of units, mean, sigma, lower 
and upper three sigma limits, minimum and maximum 
values, ete. 

At the other end of the problem scale are those 
problems which involve a relatively small input of 
data representing parameters or initial conditions, 
and a relatively complicated program to guide the 
computation of an extensive machine for a sizable 
length of time before answers are produced. Such a 
problem would be more likely solved on a large-scale 
digital computer like the IBM 704 described in the 
article by I. Silberman in this issue. Solution of such 
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complete identification will not be provided in all 
cases. 

Engineering applications of data processing ma- 
chines usually involve using more of the capability of 
the machines than those required for listing only. In 
this sense the problems range from those in which the 
information contained on one card is the basis of cal- 
culations whose results are punched on that same or 
a duplicate card. Results of equipment tests, for ex- 
ample, including voltage and current input as well as 
the useful power output, might be punched on a card. 
An IBM 604 or 650 could then be employed to caleu- 
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problems by high-speed digital computers is obviously 
not a simple task, but the results obtained justify the 
engineering effort involved. 

There are many other types of engineering prob- 
lems whose solution can contribute to scientific prog- 
ress. The broad spectrum of computing and data 
processing equipment of all types which is now, or 
will become, available provides the means of making 
equipment suitable to the scope of the problem that is 
to be solved. Even the relatively simple punched card 
equipment, by permitting the analysis of scientific 
data, can contribute to engineering capability. 
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The Bedford 
Digital Computer 
Facility 


Tue largest and fastest electronic digital computer 
now in use in Raytheon is the IBM 704 at Missile Sys- 
tems Division’s Systems Laboratory. The 704 provides 
Raytheon engineers and scientists with one of the most 
advanced systems available for the solution of mathe- 
matical problems. 

Essentially, this large high-speed digital computer 
ean add, subtract, multiply, divide, and perform logi- 
cal operations in binary form. Other operations are 
derived from those listed; for example, integration is 
performed by summation, and converging series may 
be used to approximate trigonometric functions. With 
these capabilities, any problem which ean be stated in 
logical (mathematical) terms ¢an be solved on the 
digital computer. Typical engineering problems may 
involve the solution of simultaneous algebraic equa- 
tions, the solution of ordinary or partial differential 
equations, linear programming, etc. Other desired so- 
lutions may demand the use of logic; i.., the making 
of decisions. The IBM 704 is capable of handling any 
of these problems within the limitations of its memory 
size (4096 words of core memory, 8192 of drum, and 
900,000 in each of four tape units), which is ade- 
quate for most scientific computations. 

To enable the computer to solve a problem, a de- 
tailed set of instructions called a program or a routine 
is written. The procedure of writing the program is 
one of the most time-consuming and tedious steps in 
problem solving on the computer. The problem must 
first be clearly defined, a complete analysis must be 
made of all phases of the solution, and the program 
must then be translated (encoded) into machine lan- 
guage. After these three steps are completed, the 
program, along with all data, is punched on standard 
IBM cards and fed into the computer memory units. 

The computer executes the programmed instruc- 
tions step by step, taking data from its various memo- 
ries, performing arithmetic operations, transferring 
data between the various memory units, inputing and 
outputing information, and continually executing logi- 
cal operations to do its ‘‘bookkeeping’’; i.e., keeping 
track of what has been done, how much more has to 
be done, what to do next, and when, with the program 
completed, to stop. 

The labor of writing detailed programs is consider- 
ably reduced by maintaining libraries of subroutines 
and by simplified programming techniques. The com- 
puter library at Bedford contains subroutines ranging 
from the common functions such as square root, sine, 
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cosine, logarithm, exponential, etc., to highly compli- 
cated routines for inverting matrices and solving sys- 
tems of differential equations. These routines are 
used whenever possible, requiring the programmer 
merely to adapt the engineering problem to the format 
of the library routine. The Bedford library is kept up 
to date by SHARE, the IBM 704 cooperative users’ 
association, in which the over 200 members exchange 
thousands of programs costing millions of dollars to 
write and develop. 

Another programming aid is the SAP (SHARE 
Assembly Program) assembler. Although the IBM 
704 is a binary machine, SAP allows the programmer 
to write his programs in decimal instead of binary 
form and in relative address instead of absolute ad- 
dress form; i.e., without regard for where the various 
data and instructions are to be stored in the memory. 
Mnemonic symbolic names consisting of letters and 
numbers are attached to physical quantities instead of 
identifying them by the locations where they are 
stored in the core memory. SAP assembles the pro- 
gram as written by the programmer into a binary 
program understood by the machine, converts all sym- 
bols to numbers, and properly arranges all parts of 
the program in the core memory. 

In the last few years even more powerful languages, 
the languages of automatic programming methods 
such as FORTRAN (FORmula TRANslation), have 
been developed for digital computers. In FORTRAN, 
mathematical formulae and expresssions do not have 
to be broken down into their basic arithmetic opera- 
tions, and English language statements are used to 
control the logic of the program. The FORTRAN 
compiler takes the program written in FORTRAN and 
produces both symbolic and binary programs. Al- 
though professional programmers using SAP can pro- 
duce a more efficient program than FORTRAN can 
for large complicated problems, FORTRAN is ex- 
tremely useful for getting short ‘‘one-shot’’ type 
problems done in a hurry. It is a language that can 
be learned in a few days by a person with no previous 
programming experience, and it is recommended for 
engineers who would like to do some of their own pro- 
gramming. 

Since its installation, the Bedford IBM 704 has 
been used to solve hundreds of scientific problems for 
Raytheon divisions. It is impossible to list all these 
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applications; however, some of the more important 
ones are worthy of mention. 

The majority of the work for Missile Systems Divi- 
sion has been the simulation of missile flights. The 
equations describing the trajectories of missiles are a 
system of simultaneous second order nonlinear ordi- 
nary differential equations. Many design changes in 
missile shapes, motors, and control systems have been 
arrived at by varying design parameters over the 
ranges of possible values and then computing the re- 
sulting trajectories. The optimum value of the parame- 
ter is chosen after detailed examination of the com- 
puted trajectories. 

The 704 has also been used in problems connected 
with heat flow in missiles. In the flight of a missile 
large amounts of heat are generated at the surface, 
and the rate of flow of this heat into the interior of 
the missile must be determined. This has been done 
by solution of complex partial differential equations 
relating the heat conduction, convection, and radia- 
tion, and the ablation of the surface material. 

Three Raytheon divisions, Microwave and Power 
Tube, Research, and Missile Systems, are presently co- 
operating in a project in which the 704 is being used 
to study the relationship between electron trajectories 
and electric fields in magnetrons and Amplitrons. It 
is expected that these studies will lead to improve- 
ments in design and efficiency of these tubes and re- 
veal the mechanism for the generation of noise, 
spurious oscillations, and harmonics. 

The computer technique available for the 704 which 
is growing in importance is linear programming. In 
linear programming a situation which can be ex- 
pressed by a linear function is optimized subject to 
certain restrictions. These restrictions must be ex- 
pressible as a system of linear equalities or inequali- 
ties. The number of equations describing the 
restrictions need not equal the number of variables 
to be optimized. The computer can then maximize or 
minimize the desired linear function of the variables. 
This linear programming technique is being used for 
the optimization of transistor production at Semicon- 
ductor Division and for the patching of radomes to 
reduce boresight errors for Missile Systems Division. 
Research Division also is using this technique to de- 
termine the distribution of electron emission from 
eathodes, and thus minimize the deviations from space- 
charge conditions in magnetron simulation studies. 

There are many other problems where the 704 can 
be useful and, in some cases, provide a solution not 
otherwise attainable. Raytheon engineers and scien- 
tists who may have questions on the 704 applications 
to their problems are urged to contact the author at 
the Bedford Systems Laboratory. A staff of analysts 
and programmers is available to provide assistance 
and information. 
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Solution of problem is transmitted from ana- 
lytical control unit to output units and read 
from punched cards or on-line printer. 
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Tur pace of technological development has accel- 
erated tremendously in recent times under pressures 
which we are all aware of. Because of present de- 
mands, we now have to go from the raw idea to hard- 
ware of proven reliability in only a few years or find 
ourselves in a hopeless position with respect to foreign 
and domestic competitors. We no longer have time 
for the slow evolutionary process of natural selection, 
as has been the case in the past history of technology. 

For example, one of our older precision instruments 
is the watch. A watch that is accurate to one minute 
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ferences and to vibration and acceleration in all di- 
rections. Its error is about 0.067 per cent, which ‘is 
pretty good for an instrument selling for about $10. 


Yet this is not surprising, since the first watches were 
developed before 1500. The watch of today has over 
400 years of development behind it. 

The automobile and the vacuum tube, other reason- 
ably well made mechanisms, were developed in only 
60 years. The atomic reactor, however, has been forced 
to reach a very sophisticated development in 20 years, 
as have the transistor and the various missiles. , 

This transition from idea to reliable hardware can 
be done using classical methods of development if the 
operation is staffed with geniuses in the areas of re- 
search and development, engineering, and production. 
Unfortunately, there are just not enough geniuses to 
go around. There is need for a tool which will permit 
engineers of the ordinary level of brilliance to work 
at the necessary accelerated pace. 

It is interesting that such a tool has been worked 
out by following the ordinary pattern of engineering 
work. Engineering development is always done by 
starting with a proposed model, evaluating it, and 
then making a random search for methods to 
strengthen the weak points of the design. The random 
search is needed because methods of strengthening the 
weak points are beyond the current state of the art, 
otherwise the points would have been designed 
stronger. This is the debugging phase common to all 
design projects. Ordinarily it is not recognized that 
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a random search is being conducted, since the usual 
phrasing is ‘‘We don’t know the answer yet, but we 
are working on several very promising lines of de- 
velopment.’’ The successful model with the variations 
discovered in the final phase then becomes the proto- 
type for the next model, in which design improvements 
are found by random search. 

The pattern of engineering development is very 
much like that taken in natural evolution where new 
models of a species, formed by a random combination 
of genes, are evaluated in an environment. Only those 
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steps in making a methodical search for new variants 
in well known species and evaluating them for spe- 
cial purposes. They very soon learned to interfere 
in the natural process and choose variants which they 
expected to be of interest by cross-breeding of plants 
and animals. Since natural variants of gene combi- 
nations do not ordinarily differ greatly from the 
parent stocks, they learned to interfere at the gene 
level to get new random combinations, using X-rays 
and chemicals to induce changes in the gene patterns. 
By doing this, they produced mostly monsters and 
non-viable stock, but a few extremely successful species 
emerged. (Successful is defined as being useful to 
man.) 

Having learned to produce interesting variants, the 
biologists were faced with the problem of evaluating 
them. The problem here was to find whether the new 
species would reproduce itself and survive in an en- 
vironment where it would be useful. A plant might be 
very desirable for. its fruit, but this alone could not 
qualify it as successful. It must also be able to sur- 
vive in reasonable environments. Testing of plants is 
necessarily slow because they have to be tested 
throughout a year-long growth cycle, in many climates 
and soils, and in the presence of many insect enemies 
and blights. If the plant survives at all in its original 
environment, it takes decades or centuries to get it 
thoroughly tested. So the biologists looked for a 
quicker way to get the answers. 

The break-through which made it possible to evalu- 
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ate plant species in years instead of decades came as 
a result of the work of Ronald Fisher. He developed 
methods of exposing random samples of the plant 
under test to various environments under different 
conditions of soil, rain, sunshine, fertilizer, ete., and 
of evaluating the resulting yields by statistical meth- 
ods. By using this method of evaluation, he was able 
to sort out and assess individually the factors im- 
portant to the yield. 

Fisher’s statistical evaluation was a really spectacu- 









lar achievement. Suppose a crop were planted in New 
England and another one in California. The New 
England planting produced a very poor crop (10 
units) and the California planting an excellent yield 
(25 units). Was the difference caused by the shorter 
growing season, the difference in rainfall, hours of 
sunshine, acidity of the soil, attacks by plant enemies, 
better cultivation, or what? Fisher was able to show 
how some of these questions could be answered. His 
more important contribution, however, was to point 
out that those questions which could not be answered 
caused a residual variability in yield, tending to make 
comparison of the New England yield to that of the 
California yield subject to a caleulated risk. It would 
be extremely unsound to predict that you would get 
the same yields if you repeated the experiment. It is 
well known that crop yields vary from year to year, 
and that there may be a bad yield in one part of the 
country and a good crop in another in the same year. 
Before the two yields can be compared an allowance 
has to be made for this variability. The only way to 
measure this factor is by use of statistics. 

Over a period of years, Fisher worked out methods 
of making comparisons and evaluating statistical re- 
sults. This had tremendous economic impact, and 
statistical control became a standard practice at agri- 
cultural research stations and colleges. One of the 
most interesting rules laid down by Fisher was that 
a statistician does not state that A is better than B 
unless the odds are 19 to 1 that the statement is cor- 
rect. Working with these odds, the statistician natu- 
rally makes very few wrong statements. If the odds 
are less than 19 to 1, the statistician says that he 
cannot make a statement, but that if more data are col- 
lected he may be able to make one on the basis of the 
added evidence. 

The impact of Fisher’s work was not confined to 
agricultural situations. The value for industrial con- 


trol was apparent. Much of the early work of appli- 





cation of statistical methods to manufacturing was § to 


done at Bell Telephone Laboratories. The classic text } 4 


on the field was W. A. Shewhart’s ‘‘ Economic Control quali 
of Quality of Manufactured Products,’’ published jp | this” 
1931. This text is generally held to be the first on § mp? 


quality control and laid down the definitions and fun. | t 


damental methods of operation in the field. in th 
Shewhart’s work also laid the groundwork for the qo # 
explosive growth that came in the 1940’s during the sible 





war. Under war conditions many people and compa- | idet 
nies began manufacturing products which were totally } lute 
new to them. There was no background to establish | fro 
quality levels for these products. In order to have } the 
control of quality it was necessary to define quality | ove 
and set up standards for measuring it on a nation. f thei 


wide scale. Nationwide control was accomplished by 


setting up military standards for quality measure. } whe 


ment. Since almost everything was purchased to some 
extent by the military services, quality control meth- 
ods were quickly picked up by many industries. It 
should be noted that the military standards gave very 
constant protection to the manufacturer and relatively 
weak and variable protection to the purchaser. Essen- 
tially, the manufacturer was assured that if he pro- 
dueed material of the specified quality it was almost 
certain to be accepted, but if his quality was not up 
to the standard there was still a pretty good chance 
that it would pass. These provisions were included 
largely because manufacturers became very nervous at 
the thought that they might actually have to produce 
goods at a level of quality which they were in no po- 
sition to measure. 

In order to sell to the military, the manufacturers 
had to set up quality control organizations and in- 
tegrate them into the manufacturing process, even 
though they were not forced to improve their quality 
substantially. 

Once quality control systems had been installed it 
was found that they led to great savings in inspee- 
tion and manufacturing. When the product quality 
was good, inspection costs could be reduced. When 
quality was bad, quality control methods were used 
to find the causes, which could then be corrected by 
production management. Here the savings were in 
the production cost, being the money saved by not 
making unacceptable products. 

It is important to note this dual function of qual- 
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ity control. If your product is of high quality, look 
to quality control to save inspection costs. If the 
quality of your product is less than satisfactory, use 
quality control to find the assignable causes. Given 
this information, the production force may be able to 
improve the product by eliminating the sources of 
trouble. It may be found that the trouble is inherent 
in the design, and that either design must be changed 
or the specifications altered before it would be pos- 
sible to make a satisfactory product. 

The basic ideas of quality control are very simple. 
The first principle is that it is impossible to make or 
find two objects which are identical. It is also im- 
possible to subject two pieces to exactly the same treat- 
ment. Most engineers and scientists, however, find this 
astumbling block since they believe implicitly in the 
laws of causality. They are convinced that treating 
identical pieces with the same process will lead to 
identical changes in the pieces. Logically this is abso- 
lutely correct. In practice it fails in degrees that range 
from barely detectable differences to cases in which 
the final pieces are not even recognizably alike. More- 
over, engineers and scientists have been conditioned by 
their training to look for similarities between items 
under consideration and to call the items identical 
when the differences between them become unde- 
tectable. In the most precise experiments, the engineer 
ultimately runs into the situation where his measuring 
equipment can no longer detect differences between 
the pieces. Unfortunately he often takes the philosophi- 
eal stand that anything he cannot measure does not 
exist. Since he cannot measure a difference between 
the pieces, he is likely to state there is no difference. 

This attitude might be a reasonable one to take if 
the engineer or scientist were omniscient and of un- 
limited energy. Unhappily he is seldom in a position 


to measure all the properties of the pieces on which 


he is experimenting. From the universe of properties 
which characterize each unit, he is forced to select a 
small number which he believes to be important. All 
too often, however, he denies that those he has ex- 
cluded can have any effect on the results of the ex- 
periment. 

Since he is expert in his field and careful to con- 
vince himself that his experimental design is a proper 
one, he can defend his design and show with great 
plausibility that his results are valid. Unhappily the 
Pieces being tested see all the variables which he has 
excluded from his analysis and are affected by them 
with absolutely no regard for his logical skills. 

The second basic principle of quality control is that 
if you try to make a number of pieces to a given 
measurement, most of them will come near the desired 
value and relatively few will be far from it. Under 
ordinary machine shop practice, for example, if you 
turn shafts to 1.000 in. diameter, most of them will 
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be between 0.999 and 1.001 in., a fair number between 


- 0.998 and 0.999 in. or 1.001 and 1.002 in., a few will be 


between 0.997 and 0.998 in. or 1.002 and 1.003 in., 
and almost none below 0.997 in. or above 1.003 in. 
The mean value is 1.000 in., and the range is from 
0.997 to 1.003 in. It is less likely that pieces will be 
made near the ends of the range than near the middle. 
Two concepts are required to describe the shaft: a 
measure of central tendency which tells the value 
around which the shafts tend to cluster; and a meas- 
ure of dispersion which tells how far they can be ex- 
pected to get from the central value. A group of things 
which are supposed to be the same is called a popu- 
lation. The numbers which describe the population 
are called the population parameters. One of the 
parameters describing this population is the mean, 
in this case 1.000 in. diameter. The mean is a good 
measure of central tendency. The range, however, is 
not a good measure of dispersion. As you make more 
pieces there is always a chance that you will make one 
larger or smaller than has yet been made, causing the 
range to increase, but not in uniform increments. 
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Populations of many types have been studied in 
great detail. The concept of a normal population was 
developed very early by Gauss and has been found to 
be an excellent tool in working with distributions 
which occur in practice. Suppose we take the data 
from the shaft-turning operation and plot the num- 
bers of shafts having diameters within a given in- 
erement against diameter. From 0.000 in. diameter 
to 0.996 in. diameter we will find no shafts. The 
occurrence of shafts will increase slowly as we go 
from 0.996 to 0.998 in. then more rapidly as we go 
to 0.9995 in., start decreasing as we pass 1.000 in., fall 
off more sharply at 1.0005 in., begin to decrease more 
slowly at 1.002 in., fall almost to zero at 1.003 in. and 
we will find no shafts above 1.004 in. or so. This plot 
is the characteristic bell-shaped curve mentioned in 
so many text books. The equation describing the 
curve is: 
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where x is the mean of the population, and ¢ is the 
standard deviation. (The units of o are the same as 
units of x.) This is clearly an equation which we 
would be happy not to have to work with, and as a 
matter of practice do not usually use. It can be shown 
that ¢ is a measure of dispersion of the population and 
has no connection with the mean, while x, the mean, 
has no connection with the dispersion. The mean can 
be calculated by adding up all the diameters of the 
shafts in the population and dividing by the number 


of shafts: 
1 N 
x= W > Xj. 


We can calculate ¢ by the formula: 


, N 1/2 
oi aw (x;-x) | : 
E> 


For a given set of data, X and ¢ are clearly constants 
which determine the location of the curve on the x 
axis and the distance it spreads from the mean. [f 
we subtract xX from each x value so the new X= x-x, 
this has the effect of sliding the x axis so that the 
mean of the transformed distribution is 0, but does 
not change the shape. If we now divide each trans- 
formed x value by s, giving the name t to the new X 


value (whieh is now =) , we get a new shape to the 


eurve which can now be expressed by the equation: 





This is the equation of the normalized probability 
curve. It has a mean of zero and a standard deviation 
(c) of 1, expressed in units of t. The curve is sym- 
metrical around the y axis, and asymptotic to the x 
axis on both ends. The area under this curve is unity, 
the constant 1/ / 9x having been chosen because 


1 
¥ ew dx=— yz. 
6 2 


The probability that a randomly chosen member of 
the population will be found in an interval t, to te 
is equal to the area under the curve in this interval. 
Since the total area under the curve is 1, the probabil- 
ity that any member is in some interval of t is 1. The 
probability that a randomly chosen unit lies above 
the mean (or below the mean) is 0.5. The area of the 
curve between t = 0 (the mean) and t = 1 (one sigma 
above the mean) is 0.341, and, therefore, 34.1 per cent 
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of the population lies in this interval. By symmetry 
an equal fraction lies between t = —1 (one sigma he. 
low the mean) and t = 0. Over 68 per cent of the 
population lies within one standard deviation of the 
mean. Between t = 1 and t = 2 the area is 0.136, go 
about 27 per cent lies between 1 and 2 standard devia. 
tions from the mean. The range between t = —2 and 
t = 2 eontains 95 per cent of the population. The 
range from t = —3 to t = 3 contains 99.7 per cent of 
the population, about 2.5 per cent lying in each tail 
between t = — 2 and t = — 3 andt = 2 and t = 3. 

These figures affirm our feeling about dispersion, 
that most of the population will be close to the mean 
and relatively little far removed from the mean ona 
mathematical basis. We can now tell what fraction of 
the population will be within a stated distance from 
the mean if we can find some easy way to calculate the 
mean and the standard deviation, and if we can be 
sure that the population is normal. If we look at our 
shaft-turning problem again we find that if we havea 
mean of 1.000 in. and a standard deviation of 0.001 
in., 68 per cent of the shafts lie between 0.999 and 
1.001 in., 95 per cent between 0.998 and 1.002 in., and 
only 3 in 1000 outside the limits of 0.997 in. and 
1.003 in. 

Suppose we put tolerance limits of 0.997 to 1.003 in. 
on the shaft, set up the lathe, set the mean at 1.000 
in., and start to turn. Only 3 shafts in 1000 will be 
out of limits. As turning continues, however, the tool 
starts to wear and after an hour the mean has shifted 
to 1.001 in. The standard deviation does not change, 
so our upper 2 ¢ limit is now 1.003 in., and our 
lower 3 o limit 0.998 in. We are now producing 2.5 
per cent of the shafts above tolerance and none below 
tolerance. If we let the tool wear until the mean has 
changed to 1.002 in., the upper 1 ¢ limit is 1.003 in. 
and 15 per cent of the shafts are out of tolerance on 
the high side, none on the low side. 

It is not sufficient to set the mean value of the 
shafts within the tolerance limits; we must center the | 
process so the mean is in the middle of the tolerance 
band. We cannot set tolerances rationally until we 
know the standard deviation the population will have 
when we produce it by a given process. We cannot 
produce the piece to rational tolerances unless we con- 
trol the mean at the point of production. When we 
set up to live according to these rules, we have set up 
an engineering control system and can tell what frae- 
tion defective the process should produce. If the 
fraction defective increases above the expected value, 
we can find the cause and bring the process back to 
the expected value. 

A subsequent article will demonstrate how engi- 
neering decisions can be made by statistical methods s0 
that the greatest yield of desirable product may be 
made at the lowest cost. 
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§. P. HEIMS, ‘‘Polarization of Nuclei by Macroscopic 
Rotation,’’ presented at the International Symposium on 
Polarization Phenomena of Nucleons, Basel, Switzerland, 
July 4-8, 1960. 

It has been suggested by Janes that nuclei in a crystal 
may be polarized by means of a rapid rotation of the 
erystal. This effect is small but observable. If the ro- 
tating internal magnetic field, which will act on a 
nuclear spin, is considered, it becomes clear that a 
mechanism is available for transforming part of the 
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space-charge configuration, justifying the magnetron 
hypothesis. When the power level is high, or for opera- 
tion near the edge of the frequency band, trajectories 
were found to be perturbed. 

This method of computation enabled the determination 
of the anode conversion efficiency of the Amplitron. 
Use of an analog computer to account for variation of 
electrode geometrical configuration is also discussed, 
and a means of determining the self-consistent mag- 
netron solution using adiabatic theory is presented. 


orbital angular momentum of the macroscopic rotation 
into the polarization of nuclear spins. Understanding 
of this effect is facilitated by a generalized form of 
Larmor’s theorem. The relaxation times required for 
reaching equilibrium polarization are evaluated and re- 
lated to the relaxation times found in magnetic reso- 
nance experiments involving magnetic fields instead of 
rotations, 


C. G. SHAFER, ‘‘ Design and Operation of an S-Band 
Traveling-Wave Diode Parametric Amplifier,’’ presented 
at the IRE WESCON Convention, Los Angeles, Calif., 
August 24, 1960. 

The achievement of low-noise microwave amplifica- 
tion by employing microwave junction diodes as vari- 
able capacitances has aroused interest in such devices 
by microwave engineers. These devices require a multi- 


Cc. A. KLEIN, P. P. DEBYE, ‘‘On the Statistics of Multi- 
level Systems in Semiconductors,’’ presented at the Inter- 
national Conference on Semiconductor Physics, Prague, 
Czech., August 29, 1960. 


frequency cavity or waveguide system in which the 
signal frequencies are coupled by a time-varying ¢a- 
pacitance. Since the product of the voltage gain and 
fractional bandwidth is a constant for a given system, 


The electron population statistics of the multilevel 
system associated with a substitutional gold impurity 
in germanium have been recently investigated by 
Khartaiev. He computed the probability, P, of having 
from zero to an infinite number of electrons trapped by 
the gold ‘‘core.’’ Expressions similar to Khartaiev’s 
have been used by Mashovets, in the case of copper-doped 
material. We have examined this question from the 
grand canonical ensemble point of view, and have shown 
that the Khartaiev formulas cannot be sustantiated by 
a rigorous treatment. Our computations for the case of 
p-type germanium were in agreement with the result of 
a straightforward application of the law of mass action. 
Moreover, the approach used provides a correct interpre- 
tation of the statistical weight factors. This formula- 
tion has been applied to carefully evaluated low-temper- 
ature carrier concentration data with the aim of arriv- 


achievement of high gain is accomplished at the loss 
of bandwidth. A method of overcoming this restriction 
is to conceive a structure which has a large bandwidth 
and low gain at each stage, but cascades stages to 
achieve a large over-all gain and bandwidth. Subject 
to the design requirements and by use of the ‘‘ coupled- 
mode’’ theory, an amplifier was designed, which has the 
following performance characteristics: 


Frequency range 3240-3400 me 
Forward gain per diode ~1.0 db 
Pump power (for 4 diodes) 400 mw 
Pump frequency 6600 me 
Reverse loss per diode ~ 0.25 db 


M. G. HOLLAND, ‘‘Low Temperature Thermal Conduc- 
tivity of Silicon,’’ presented at the International Confer- 
ence on Semiconductor Physics, Prague, Czech., August 
29, 1960. 





ing at an improved description of the donor level and 
of the first acceptor level. Evidence has been accumu- 
lated which demonstrates that Hall effect experiments 
indeed allow a rather accurate estimation of the statisti- 
cal weight factors, 


G. E. DOMBROWSKI, ‘‘Large-Signal Analysis of the 
Amplitron,’’ presented at the AIEE-IRE Tube Research 
Conference, Seattle, Wash., June 29, 1960. 
This paper presents some experimental evidence sup- 
porting the theory that the space-charge behavior in an 
Amplitron operating at low or moderate power levels is 
essentially that of the magnetron. Large-signal trajec- 
tory calculations were made on a digital computer using 
the space-charge model for the space-charge potential 
and a large number of Laplacian rotation waves for the 
Amplitron ‘‘cireuit potential.’’ The trajectories cal- 
culated for usual Amplitron operating conditions showed 
that the non-synchronous waves have little effect on the 


The low temperature thermal conductivity of large 
samples of single-crystal silicon has been measured 
between 4° and 300° K. Samples grown by the floating- 
zone process were found to have a maximum thermal 
conductivity two times greater than samples grown from 
a erucible. As the boron content was increased from 
1013 to 5x 105 per ec, there was a 20 per cent increase 
in the maximum conductivity of the floating-zone mate- 
rial. This was somewhat dependent on the presence of 
compensating group V impurity atoms. In the crucible- 
grown material the maximum thermal conductivity was 
found to be further increased by high temperature an- 
nealing at 100° C; however, the increase was related to 
the dislocation density in the material. The effect of 
annealing at lower temperatures is also presented. 
Some correlation of the maximum in thermal conduc- 
tivity and dissolved oxygen content is shown, and evi- 
dence which identifies oxygen as the scattering center is 
discussed. 
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New Raytheon Cathode Ray 
Printer Tubes Print as Fast 
as Computers Think! 


Raytheon cathode ray printer tubes now make possible the design of extremely 
high speed electronic printing equipment capable of printing 20,000 characters 
per second — or over 10,000 lines of computer information per minute — fast 
enough to keep pace ss the ee of modern computers. 


These new gubes ete into fotuy: ‘Ahost of amazing new applications. Instant 
coast-to-coast tranismisdlog of mail grontinuous photographic recording of 
freight traffic movement — “incredibly “fast printing of labels, or facsimile data 


— these are just a few examples. 2 


Four developmental Raytheon printce tubes, the CK 1366 and CK 1368 with 
single row printing elements and the forty row CK 1367 and CK 1369, are 
available for producing writing widths of approximately 2%4 and 8% inches. 
For complete technical data and design assistance, please write to: Raytheon, 
Industrial Components Division, 55 Chapel St., Newton 58, Mass. 


RAYTHEON COMPANY 


INDUSTRIAL COMPONENTS DIVISION 


HOW IT WORKS. Raytheon printer tubes 
employ a high-resolution, high-current 
gun which sweeps a printing element con- 
sisting of a row of 250 wires per inch. 
Electron charges in the form of the de- 
sired characters are produced by deflec- 
tion and amplitude modulation of the 
beam. The charges are established on the 
paper which acts as a dielectric as it 
passes before the faceplate of the tube. 
The characters are made visible through 
the use of a positively charged ink and 
a fixing process. Since the electron 
charge on the paper is proportional to the 
beam current, fine halftones can also be 
obtained. 
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